Global Phase Space of Coherence and Entanglement in a Double- Well BEC 
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Ultracold atoms provide an ideal system for the realization of quantum technologies, but also for 
the study of fundamental physical questions such as the emergence of decoherence and classicality 
in quantum many-body systems. Here, we study the global structure of the quantum dynamics of 
bosonic atoms in a double- well trap and analyze the conditions for the generation of many-particle 
entanglement and spin squeezing which have important applications in quantum metrology. We 
show how the quantum dynamics is determined by the phase space structure of the associated 
mean-field system and where true quantum features arise beyond this 'classical' approximation. 
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Maintaining and controlling the coherence of quan- 
tum systems over time is one of the major challenges in 
contemporary physics. Low-temperature quantum gases 
trapped in optical lattices are an important instance of 
this challenge, for they provide versatile testbeds both 
for idealized models of exotic solid state systems and 
for applications in quantum optics and quantum infor- 
mation processing [IHl]. A Bose-Einstein condensate 
(BEC) loaded in two linearly coupled wells, called a 'BEC 
dimer' or 'Bosonic Josephson Junction' [5j is a partic- 
ularly appealing case, as its quantum behavior is both 
theoretically tractable [6H9] and experimentally accessi- 
ble. In particular, BEC dimers have been shown to al- 
low coherent manipulation of quantum states even on 
an atom chip, enabling matter- wave inter ferometry [10] , 
long phase coherence times and number squeezing [TTT — 
[13], and proposed chip-based gravity detectors [14J. Re- 
lated questions have been addressed in the context of the 
mathematically equivalent Lipkin-Meshkov-Glick model 
[15J- 

For large atom numbers, the coarse dynamics of the 
BEC dimer is excellently described by a 'classical' mean- 
field approximation. Recent experiments have precisely 
mapped out the emerging classical phase space struc- 
ture fT6l [T7] , showing macroscopic quantum oscillations 
as well as the emergence of self-trapped states through a 
'classical' bifurcation. Recent theoretical studies have fo- 
cused on strategies to prepare highly entangled Einstein- 
Podolsky-Rosen (EPR) states in the BEC dimer [TM2T] . 
These truly quantum phenomena appear to depend on 
the structure of the associated classical phase space [6- 
[9j [18] , but what is the exact correspondence? 

Here, we introduce a global phase space (GPS) picture 
of the dynamics of the quantum dimer, focussing on en- 
tanglement, spin squeezing and decoherence of the atoms. 
Our global picture shows that both the time evolution of 
entanglement and coherence depend very fundamentally 
on the initial coherent state. In particular, we find that 
self-trapping (ST) excellently supports many-particle en- 
tanglement. To link quantum and classical aspects of the 
dynamics, we approach the problem from three perspec- 



tives. First, for the full quantum problem, we employ 
the two-site Bose-Hubbard model and solve for all of the 
above quantum observables as functions of time. Second, 
we recall that the standard mean-field treatment of the 
dimer - the two-site Gross-Pitaevski equation - corre- 
sponds to an integrable classical dynamical system, and 
we discuss the nature of the trajectories in this system. 
Third, we use a recently developed semiclassical 'Liou- 
ville dynamics approach' (LiDA) [22J, which we argue 
is intermediate between the fully quantum and the clas- 
sical systems. Comparing these three perspectives, our 
GPS approach provides an intuitive understanding of the 
quantum observables mentioned above and insights into 
true quantum effects beyond the mean-field or semiclas- 
sical approaches. 

Quantum and semiclassical dynamics.- The coherent 
dynamics of the Bose-Einstein condensate dimer is given 
by the two-mode Bose-Hubbard Hamiltonian (BHH) 

H = -J(a\a 2 + a^ai) + y (a^a\ + a^a^j (1) 

where dj denotes the annihilation of one bosonic atom in 
state j. The modes can be either realized by two sites 
in a double- well trap [23] 24J) or two internal states of 
the atoms [17] . The tunneling rate J and the on-site 
interaction U can be tuned individually, e.g., via a Fes- 
hbach resonance or by changing the depth of atom trap 
[U G3J [T7J [25H27] . Throughout this paper we set h = 1, 
thus measuring all energies in frequency units. If not 
stated otherwise, we assume that J = 10s _1 , which is a 
common setting in ongoing experiments [13J H7J [27] , and 
we assume a total number of N = 40 atoms. 

The mean- field counterpart of the Bose-Hubbard dimer 
is a two-mode Gross-Pitaevskii or nonlinear Schrodinger 
equation. This equation can be rewritten as an integrable 
classical Hamiltonian system with the Hamiltonian [28] 

Az 2 / 

H d (z, cp) = — V 1 - z 2 cos((p) . (2) 

The mean-field (classical) phase space consists of the two 



1 2 1 2 '0 1 2 "°- 5 

<p[tt] (p[ir] (p[n] 

Figure 1: Global phase space structure of the Bose-Hubbard dimer for (a) A = 5, (b) A = 1.5, (c) A = 0.5. The color code 
shows the condensate fraction c T (ip,z) at r = 1 s as a function of the initial state \<p,z). Lines of equal condensate fraction 
mimic the classical trajectories (black lines), despite the clear qualitative differences for the three values of A shown. However, 
we observe strong deviations in the coherence properties of the quantum state, such as S-symmetry breaking and enhanced 
condensate fraction near the unstable fixed point Fi — (tt, 0), which are clearly beyond the classical mean- field description. 



conjugate variables, which following previous usage [28] , 
we have chosen to be the relative phase ip = ipi — cp 2 £ 
[0, 2tt) between the two wells and the population imbal- 
ance z = (Ni — N 2 )/N £ [0, 1], where TVi and N 2 denote 
the number of atoms in each well and N = N± + N 2 is 
the total atomic population. The classical trajectories 
simply follow the lines of constant (conserved) energy 
Hd(z,(p) = const. They are determined by the initial 
position (po,zo) in phase space and the ratio of the in- 
teraction and the tunneling energy A = UN/ (2 J). 

The Gross-Pitaevskii equation implicitly assumes a 
pure BEC at all times. Some quantum features be- 
yond this rough approximation, in particular the quan- 
tum mechanical spreading over time, can be included if a 
quantum state is represented by a (quantum) phase- 
space density such as the Husimi function Q(jp,z) = 
\(cp, z\(po, ^o)| 2 instead of a single trajectory. Here, \ip,z) 
denotes an atomic coherent state [55], which is noth- 
ing but a pure BEC. As shown in [22J, the dynamics of 
the Husimi function follows a classical Liouville equation 
with the Hamiltonian Q plus quantum correction terms 
vanishing as 1/N. In the semiclassical LiDA [22J we thus 
represent an initial state by an ensemble of trajectories, 
whose initial positions are distributed according to the 
Husimi function of the initial quantum state. 

Global phase space structure.- We shall analyze the 
global phase space structure of the Bose-Hubbard dimer 
with special respect to its classical and quantum prop- 
erties. Therefore, we consider the dynamics of an ini- 
tially pure BEC as a function of the parameters zq and 
(fo with a focus on the purity of the condensate and the 
emergence of entanglement and spin-squeezing. The pu- 
rity is measured by the condensate fraction c T defined as 
the maximum eigenvalue of the reduced single-particle 
density matrix (SPDM) Pij = {a\dj)/N [9l [22J [30J . A 
related GPS approach was introduced in terms of the 
phase space entropy [31 J, which provides great insight 
into global dynamical properties of the system. 

The global phase space structure of the Bose-Hubbard 



dimer is shown in Fig. [T] where the condensate fraction 
c T ((p,z) at time r = Is is plotted as a function of the 
initial state \<po,zo) f° r three different values of A. The 
corresponding classical mean-field dynamics is overlaid 
as solid black lines. For A < 1, the atoms show sim- 
ple Rabi oscillations between the wells with stable fixed 
points at F\ = (0,0) and F 2 = (tt,0). According to 
the average phase (p, the oscillations are referred to as 
'zero-phase' or plasma oscillations ((p = 0) around F\ and 
'7r-phase' oscillations {(p = tt) around F 2 . The classical 
dynamics undergoes a bifurcation at A = 1, separating 
the Rabi (0 < A < 1) and Josephson (A > 1) regimes. 
The fixed point F 2 becomes hyperbolically unstable for 
A> 1, bifurcating into two self-trapping (ST) fixed points 

at Fst = (tt, zst), where zst = ±yjl - 1/A 2 |m In ad- 
dition to 7r-phase ST with average phase tt (Fig. jib), also 
'running phase' ST appears for A > 2 (see Fig. Uji). 

Comparing the condensate fraction c T (ip,z) with the 
mean- field dynamics we find that lines of equal c T ((p,z) 
mimic the classical trajectories in many respects. The 
GPS mapping for c clearly reflects zero- and 7r-phase os- 
cillations and the stable fixed points (Fig. [ljc)) as well 
as ST (Fig. [TJa-b)). A drastic drop in the condensate 
fraction is observed in particular near the classical sepa- 
ratrices, indicating a deviation from a pure BEC and thus 
the failure of a simple mean-field description. Notably, 
this loss of coherence is not strongest at the unstable 
fixed point, cf. Fig. [TJa-b). Furthermore, the quantum 
dynamics breaks the symmetry S : (ip,z) — >• ((/?, —z), as 
Cr(S(cp, z)) 7^ c T (ip,z). In the classical limit S is pre- 
served: H(S(<p,z)) = H((p,z), see Eq. ^. Of course, 
both the classical and the quantum dynamics are sym- 
metric with respect to (<p,z) — >• (—ip,—z), which trans- 
lates to relabeling wells 1 <H> 2 without changing the ini- 
tial state. 

The correspondence of quantum and classical phase 
space structure is more than a qualitative coincidence. In 
Fig.[2](a,b) we compare the quantum results (based upon 
the BHH) for the minimum condensate fraction with the 
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Figure 2: The GPS approach reveals novel features of quan- 
tum dynamical origin beyond the classical mean-field descrip- 
tion and the semiclassical LiDA. (a-b) Sections through the 
GPS for A = 5 comparing quantum results (based upon the 
BHH, red curves) with LiDA (green curves). The mean- field 
(MF) description assumes a pure BEC (blue dashed line), 
(a) Shown is the condensate fraction min(c T (0.857r, z)) in the 
time interval < r < 0.5 s. The point (0.857r,0) is near the 
hyperbolic fixed point. The LiDA reproduces ^-symmetry 
breaking, but can distinctly deviate from the quantum results 
near the hyperbolic fixed point and the separatrices (dashed 
vertical lines), (b) For min(c T (0, z)) the LiDA deviates from 
the quantum dynamics for \z\ < 0.8, i.e. near the separatri- 
ces (which are located at z — ±0.8, dashed lines) (c) Large 
condensate fraction c T in the vicinity of the classical fixed 
points is reproduced qualitatively by a time-independent mea- 
sure. The color code shows the maximum overlap of the initial 
state with a stationary state, A((p, z) — max n \ ((p, z \ E n )\ for 
A = 1.5. However, several quantum features depicted in the 
GPS, such as the symmetry breaking, are not reproduced by 
this measure, indicating that these are of dynamical origin. 



prediction of the semiclassical LiDA. The good agree- 
ment of both quantities reveals that the loss of quantum 
coherence c r (<£, z) can be mostly attributed to the classi- 
cal spreading of Husimi function. Around a stable fixed 
point neighboring trajectories remain close for all times, 
such that there is no spreading. 

But why does the quantum coherence remain reason- 
ably high near the unstable fixed point? We recall that 
Eq. (§ can be mapped to a pendulum with variable 
length l(z) = yl — z 2 and angular velocity z = ip. 
Within that picture, the unstable fixed point F 2 = (ir, 0) 
is reached when the pendulum is in its upright position 
with no angular velocity. Close to the unstable fixed 
point, the classical dynamics slows down asymptotically 
("freezes"); therefore spreading near the unstable fixed 
point is slow. Moreover, when the trajectory is mov- 
ing towards F 2 the dynamics slows down, while moving 
away from F 2 the pendulum accelerates, which breaks 



the ^-symmetry. On the contrary, close to the sepa- 
ratrix, the classical dynamics leads to derealization in 
phase space which induces decoherence of a many parti- 
cle state [221 151]. 

To study to what extent these features can be cap- 
tured in a time-independent framework, we analyze how 
much a state \tp,z) overlaps with an eigenstate \E n ) of the 
Bose-Hubbard dimer. Figure [2] (c) shows the maximum 
overlap A max ((^, z) — max n z | E n )\ for the same pa- 
rameters as in Fig. [I] (b). The overlap with the eigen- 
functions is high at the stable fixed points and minimal 
along most parts along the separatrix. This reveals one 
possible mechanism to guarantee long-time coherence of 
a BEC: If the overlap approaches unity, the initial state 
is almost stationary such that the condensate fraction re- 
mains virtually constant. This is consistent with (but the 
converse of) the results reported in refs. [22, 31J, where 
it was shown that derealization in phase space induces 
decoherence of a quantum state. Moreover, there is as 
well an eigenfunction localized around the unstable fixed 
point F 2 , which explains the surprisingly slow decoher- 
ence of the quantum state localized at this point. 

Entanglement.- A distinguishing feature of experi- 
ments with two-mode BECs is that the quantum state 
of the atoms can be manipulated with astonishing pre- 
cision. In particular, the atoms can be strongly entan- 
gled, which enables important applications in precision 
quantum metrology [T3J EH 133] • For the important spe- 
cial case of EPR entanglement, a simple criterion reads 
E > 0, where E = \(a\d 2 }\ 2 - (ajaia^) US HO]. Pro- 
ceeding in the same way as for the condensate fraction, we 
obtain a GPS picture for the EPR entanglement shown 
m Fig.[3](a,b). The measure E T ((p,z) again excellently 
mimics the classical phase space trajectories, however, 
an actual entangled state (i.e., E > 0) is found at time 
r = Is only near the classical ST fixed points (71", z$t)- 
Movies showing the temporal evolution of E T (ip,z) for 
r G [0, 3] s can be found in the supplementary material. 
Strikingly, the GPS for entanglement measured at differ- 
ent times r unveils the following behavior: If entangle- 
ment appears on timescales long compared to the period 
of the plasma oscillations, then entanglement is concen- 
trated around one or more fixed points in phase space, 
which, surprisingly, includes the unstable fixed point F 2 . 

Given the apparent similarity of the GPS for E and 
the condensate fraction c, we ask, how is E related to c? 
By approximating (ajaia^) ~ (a\a\) {a\a 2 ) we find a 
semiclassical measure E sc = |(aja2)| 2 — (a\ai) (a 2 a 2 ) = 
N 2 (c 2 - c). In Fig. [^c) a GPS for E sc is shown for the 
same parameters, with remarkable qualitative agreement 
with the exact quantum results for E. In Fig. (3^d) we 
report the temporal evolution of E and E sc of the state 
|0.2,0) which shows surprisingly good agreement except 
for an offset. However, as E sc < by definition, it does 
not serve as an entanglement criterion. 

Spin squeezing indicates a form of many-particle en- 
tanglement which is particularly important for quan- 
tum metrology fT3] 132]. A state is spin-squeezed if the 
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Figure 3: Global phase space structure of EPR entanglement 
and spin squeezing for A = 5. (a) EPR entanglement E > 
is found at r = Is solely around the ST fixed points. (b,c) 
The entanglement parameter E T (tp,z) (b) and the semiclas- 
sical approximation E sc = N 2 (c 2 — c) (c) at time r = 1 s are 
shown as a function of the initial state \<p,z). While these 
snapshots depict the GPS at r = Is, the movie in the sup- 
plementary material shows the time evolution of the GPS in 
timesteps of 0.025 s. (d,e) Time evolution of EPR entangle- 
ment (d) and spin squeezing (e) for an initially pure BEC 
\<p,z) = |0, 0.2). The semiclassical approximation reproduces 
the overall behavior quite well but exhibits an almost con- 
stant shift. A revival of entanglement is observed at t ~ 3 s. 
Spin squeezing (£ 2 < 1) is found to be in excellent agreement 
with the EPR entanglement result (E > 0) at t « 3 s. 

quantum uncertainty in a Bloch sphere representation 
is smaller than that of an atomic coherent state, i.e., 
a pure BEC. This representation is defined via the op- 
erators J x = ^(a^ai + aj<22), J y = ^(a\a\ — a\a2) , and 
J z = ^(d\ai — ^2^2), which form an angular momentum 
algebra with quantum number J = N/2 [22j [29j [30J . A 
quantum state is spectroscopically squeezed if the param- 
eter [33J 

e := N{ ^ ? (3) 

is smaller than one. Here, J n . =rii ■ J is the projection of 
the total angular momentum operator J onto n^, where 
^1,2,3 are mutually orthogonal unit vectors and AJ n is 
the variance of J n . In Fig.pne) we report the dynamical 



evolution of £ 2 for an initially pure BEC 1 0.2, 0) in excel- 
lent agreement with E(t) (d), including the entanglement 
revival at t ~ 3 s. 

Discussion.- While chaotic classical dynamics has been 
shown typically to lead to a fast (exponential) decay 
of the coherence of a many-particle state [34J, little 
is known about the coherence of quantum self-trapped 
states (especially in a high dimensional phase space), 
where the classical counterparts are (quasi-) periodic or- 
bits, although some insight can be provided by semi- 
classical approaches beyond the mean- field limit [35H37] . 
Fig. [2] reveals one aspect of the coherence of self-trapped 
states. For z near ±1 the many-particle state is close 
to a pure condensate at time r, whereas the conden- 
sate fraction drops drastically well before the separatrix 
(dashed line) is reached. In Fig. [lja-b) the dark red re- 
gions around the fixed points Fst further confirm that 
the self-trapped regions maintain high condensate frac- 
tion. Hence, ST is a mechanism to preserve coherence of 
many-particle states. In contrast, in the vicinity of the 
separatrix the condensate fraction falls off sharply. 

We have analyzed the connection between quantum ob- 
servables of the BEC dimer and the structure of the un- 
derlying classical phase space, including fixed points, sep- 
aratrices, and ST. The question remains: how well do our 
results carry over beyond the Bose-Hubbard model? Re- 
cent numerical studies beyond the Bose-Hubbard model 
(including higher-lying states in the individual wells) 
[SHI EH] show that ST is only present as long as the sys- 
tem remains coherent. Hence, there is evidence that our 
results do reflect a fundamental relation between ST and 
coherence of Bosonic quantum systems. 

In larger optical lattices, the self-trapped states of the 
dimer correspond to 'discrete breathers' or 'intrinsic lo- 
calized modes' [40-42J). These correspond to classical 
trajectories which are practically embedded on a two- 
dimensional torus in the high dimensional phase space 
and are thus (quasi-)periodic in time [4TH44] . In other 
words, discrete breathers involve localization in phase 
space. Moreover, discrete breathers become attractive 
fixed points in presence of dissipation and have been 
found in a variety of physical and biochemical systems 
[40| [42] . Just as derealization in classical phase space 
leads to decoherence [22j [31] . we expect that discrete 
breathers are candidates to support long-lived coherent 
and possibly entangled many-body states even in com- 
plex dissipative systems, such as biomolecular systems. 
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